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^Sl , Abstract 

^~>. How far is neuroepithelial cell proliferation in the developing central nervous sys- 

tem a deterministic process? Or, to put it in a more precise way, how accurately 
can it be described by a deterministic mathematical model? To provide tracks to 
answer this question, a deterministic system of transport and diffusion partial dif- 
ferential equations, both physiologically and spatially structured, is introduced as a 
model to describe the spatially organized process of cell proliferation during the de- 
velopment of the central nervous system. As an initial step towards dealing with the 
three-dimensional case, a unidimensional version of the model is presented. Numerical 
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r^ , analysis and numerical tests are performed. In this work we also achieve a first ex- 
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perimental validation of the proposed model, by using cell proliferation data recorded 
from histological sections obtained during the development of the optic tectum in the 
chick embryo. 



> ■ 1 Introduction 

^^ ' The present work aims at modeling the spatial organization of the neuroepithelial 

(NE) cell proliferation in a developing cortical structure along an early and brief 
developmental period by using sets of quantitative data empirically obtained from a 
standardized experimental model: the developing chick optic tectum (OT). The math- 
ematical model is based on a deterministic approach that uses the formalism of partial 
differential equations. 



rS I 1.1 Biological background: Developmental neurobiology 



1.1.1 Relevance of a spatial and temporal organization in a devel- 
oping system 

The appropriate number of cells of each terminally differentiated cell type and also 
the spatial patterns they exhibit within the different tissues and organs composing 
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pluricellular organisms are governed by interactive self-regulating behaviors that the 
developing cells exhibit during the embryonic development Jj . The increase in supra- 
cellular complexity generated during development requires the temporally and spa- 
tially organized operation of specific developmental cell behaviors (DCBs). These 
DCBs are usually reciprocally regulated and operate simultaneously and interactively 
[5]. Every developing cell population can be considered as both emitter and receiver 
of developmental regulatory signals having, on the one hand, informative and, on the 
other hand, structural roles. Thus, a central hypothesis in Developmental Biology 
proposes that the space-time organized operation of developmental cell behaviors de- 
pends on the cooperative establishment of spatially organized cell signaling networks 
mediated by diffusing informative molecules [31 SlEl E] ■ Molecular diffusion results in 
asymmetric distribution of developmentally active informative signals. This asymme- 
try plays a fundamental role in establishing temporal and/or spatial organization of 
specific DCBs that result in the whole developmental process in the organized pat- 
terns of cells, tissues and organs, that living organisms eventually exhibit in their final, 
terminally differentiated, state. 



1.1.2 The organizers and CNS patterning 

Not all developing cell populations possess similarly relevant informative roles. There 
exist specific transient cell populations, the so-called organizers, that primarily play in- 
formative roles influencing or regulating the developmental behavior of the other cells. 
By means of installing asymmetric distributions of developmentally active signals the 
organizers serve as instruments of an informative reference system in the establish- 
ment of spatially organized processes of cell determination and differentiation [71 [S] . 
Amongst the most complex biological structures, the multilayered concentric neuronal 
organization of the central nervous system (CNS), i.e., brain cortex, cerebellum cor- 
tex etc., occupy a privileged position. The development of such a structural and 
functional complexity, the so-called corticogenesis, requires the organized operation of 
several DCBs. Amongst these DCBs: (a) the cell proliferation (CP), counteracted by 
apoptosis or programmed cell death, is involved in the generation of the appropriate 
number of neurons for each cortical area and each cortical layer; (b) the directed cell 
migration, a process mediated by specific interfacial interactions between cell surface 
and extracellular matrix components, controls the correct position of each specific 
neuronal type along the CNS spatial axes; (c) cell determination or commitment, a 
process mediated by irreversible genetic information reprogramming, allows different 
cell populations to select one out of a set of multiple developmental pathways; (d) cell 
differentiation, a process mediated by selective gene activation and selective protein 
synthesis, warrants the expression of different specific neuronal phenotypes; finally, 
(e) the development of the neural processes (neuritogenesis) , dendrites and axons and 
the establishment of specific synaptic contacts (synaptogenesis) conclude this complex 
self-organizing and interactively regulated process of corticogenesis. 

All these DCBs are properly organized in time and space. During the early stages, 
several organizers strategically positioned along the cephalic-caudal and dorsal- ventral 
axes of the developing CNS establish the primitive pattern of each region and sub- 
region, as well as the identity of the neuromeres. Later, cell proliferation, migration 
and differentiation controlled by reciprocal interactions allow the whole system to ex- 
pand and differentiate into several structurally and functionally integrated systems of 
interconnected neuronal circuits. Abundant evidence indicate that all these processes 
are spatially and temporally organized [SI [TUl El [TH [T31 [T3] . 



1.1.3 Roles of neuroepithelial cell proliferation and neuronal mi- 
gration in corticogenesis 

During the early development, the primitive CNS primordium, the so-called neu- 
ral tube (NT), is almost exclusively composed of proliferative neuroepithelial (NE) 
cells. During this early proliferative phase the NE cells behave as a population of 
self-renewing stem cells that divide symmetrically, i.e. from each dividing NE cell 
originates two similar NE cells. These dividing NE cells are typically located along 
the inner surface of the NT forming the so-called generation zone (GZ). This type 
of early proliferation produces an amplification of the NE cells: they significantly in- 
crease in number, and the GZ increases in both area (planar expansion) and thickness 
(radial expansion). During a later proliferative phase NE cells proliferate asymmetri- 
cally, i.e. from each dividing NE cell originates a NE cell and a postmitotic neuron 
(PN). The number of NE cells and the GZ size stabilize while postmitotic neurons ac- 
cumulate at an underlying premigratory zone (PMZ). Later, neurons migrate radially 
and the PMZ weakens while different layers of concentrically organized neurons begin 
to organize the future cortex. By the end of the proliferative phase the number of NE 
cells decreases, the GZ becomes thinner and finally disappears. 



1.1.4 Neuroepithelial cell proliferation analyzed as a stochastic point 
process 

The NE cells compose a synchronized and functionally integrated population of stem 
cells. Biomathematical analyses of numerical series representing the position of mitotic 
NE cells along the OT cephalic-caudal axis, analyzed within the framework provided 
by the theory of stochastic point processes, show that NE cells do not proliferate in- 
dependently. It can be shown that the proliferative dynamics embeds (subsumes) two 
different components: (a) a non-stationary component that can be deterministically 
described and reveals a spatially organized process and (b) a stationary stochastic 
point process representing a uniformly distributed basal proliferative activity I15[ 116) . 
More recently a different approach allows demonstrating that the stochastic compo- 
nent corresponds to an anti-correlated stationary process implying the result that 
there are short-range inhibitory interactions between neighboring proliferating NE 
cells [17j [18] . These mathematical results are consistent with experimental evidence 
indicating that cell proliferation is a DCB submitted to cell-cell controlling interactions 
and also to intracellular controlling processes. 

1.1.5 Extra- and intracellular signals controlling proliferation dy- 
namics 

A general notion establishes that the coordinated behavior shown by proliferating cells 
during development |19) depends, on one hand, on a set of extracellular signals, the 
so-called growth factors [20j [21] and a set of intracellular informative signals compos- 
ing a cell cycle control center whose main members are the regulatory proteins cyclins 
as well as cyclin-dependent kinases [^ . 

The temporally organized operation of these regulatory proteins installs the so- 
called cell division cycle (proliferative cycle). The cell division cycle is composed of 
a periodic ordered sequence of cyclic intracellular molecular events with typical du- 
rations. It involves two main phases: (a) the mitotic or M phase that includes the 
karyokinesis (nuclear division) and the cytokinesis (cytoplasm division) and (b) the 
interphase, or intermitotic phase. In turn, each of these phases involves temporally 
ordered sub-phases. The M phase is composed of the pro-, prometa-, meta-, ana- and 
telo-phase while the interphase is composed of the Gl, S and G2 phases. Each type 
of proliferative cells population displays typical cell cycle length oscillating around a 
mean value. The total cell cycle length and also those of the different phases, vary. 



depending on the cell type. Besides, particular cell types may change their cell cycle 
durations when long time series of subsequent periods are considered. 



1.2 Biological settings of the present study 

1.2.1 The experimental model system: NE cell proliferation in the 
developing avian optic tectum 

The avian optic tectum (OT) is a mesencephalic alar plate derivative. The main OT 
input proceeds from the retinal ganglion cell axons. This information is integrated 
into others inputs, processed and transferred to other CNS areas. The OT possesses 
a typical cortical, multilayered concentric organization, consisting of alternating neu- 
ronal and fibrous layers. The developing avian OT is a well-characterized experimental 
model of CNS development. It is known that the OT determination and the positional 
information specifying the OT patterning are submitted to the organizing influence of 
the Isthmic Organizer (IsO). The IsO organizing activity is mediated by informative 
molecules which specify positional information along the OT cephalic-caudal axis [9, 
10, 11, 12]. The positional information is firstly established as asymmetric gene ex- 
pression and is then "translated" , by the differential operation of several DCBs into a 
gradient of cytoarchitectonic differentiation. There are well-documented data about 
its proliferative kinetics [221 HSl [HI [H] , the postmitotic neuronal migration [13] and 
the expression of enzymes involved in neuronal migration [11] , the histogenesis of its 
multilaminated architectural pattern |10| and also on the temporal and spatial de- 
velopmental pattern of its specific neuronal types [131 121]. There are also detailed 
structural, pharmacological and molecular studies on the development of its extrinsic 
innervation (serotonin innervation, retino-tectal connections) [S] 1251 126) . The devel- 
opmental patterns of expression of several molecular components involved in synap- 
togenesis of the OT local circuits are also well characterized [TH [371 |5S] . 



1.2.2 Biological methods and cell proliferation records 

Empirical data on NE cell proliferation analyzed in this paper were recorded from com- 
plete histological serial sections obtained from OT of different developmental stages. 
OTs obtained from 2, 4 and 6 embryonic (E) days {E2, Ei, E6) were used. During 
this interval the OT undergoes significant changes in size and shape. After dissection, 
OTs were processed for conventional histological methods [10]. Specimens were spa- 
tially oriented in order to obtain planar sections fulfilling criteria of adequacy between 
planar section orientation and developmental gradient axis position. These criteria 
were given in [6] . Cell proliferation records (CPR) along the OT cephalic-caudal axis 
were constructed on a video camera screen (Axioplan 2 imaging optical epifluorescence 
microscope coupled to an Axiocam HR color digital scanner and a computer equipped 
with the software Axiovision (all Carl Zeiss, Germany)) at a final magnification of 
1000 X. A CPR is defined as a numerical sequence that indicates the density of mi- 
totic cells in successive ordered, 25/im length windows, along OT cephalic-caudal axis. 
The total number of cells in the GZ was also computed for each 25/^r7i window. 



2 Complete mathematical model 

The present modeling concerns the development of the OT between embryonic days 

E2 and E6. It only describes the elongation of the cephalic-caudal axis and thus is 

only one space dimension. It also assumes that: 

(a) the NE cells divide symmetrically (the number of asymmetric division is negligi- 

ble), 



(b) the number of postmitotic neurons in the GZ is neghgible and the increase in cell 
number in the GZ is mainly due to the symmetric divisions and 

(c) the cell cycle length does not change significantly neither along the interval £'2- 
E4: nor between EA and E6. Consequently the numerical analysis and experimental 
validation treat separately these two periods. 

We follow a classical description of the cell cycle that can be found in f29: and 
references therein. Our model is based on equations for the cell population densities 
structured by their progression along each phase < a < 1 and their position x along 
the cephalic-caudal axis. The cells are divided in subclasses, as follows: 

ni{t, a, x) =linear density of neuroepithelial cells in phases Gi, S, G2, 

712 (i, a, x) =linear density of neuroepithelial cells in phase M, 

713 (t, a, x) = postmitotic neurons generated by asymmetric mitosis. These cells stop 
proliferating; they are destined to migrate and later differentiate into specialized neu- 



We propose a unique model for time t between days E2 and £^10. Cells in the 
phases are also subject to random motion along the cephalic-caudal axis and this 
leads us to represent the dynamics of cell division and motion by the equations 

aril I d[vi{t)ni] _ a^[D(x)m] 
at "•" da ~ dx^ ' 

8«2 I a[f2(t)»2] _ d'^[D(x)n2] (l) 

at "*" da ~ Ih? ' ^ ' 

^ ^ q{t,x)v2{t)n2it,a ^ l,x). 

This system is completed with boundary and initial conditions 

V2{t)n2{t, a — 0,x) ^ vi{t)ni{t, a — 1, x), 

Vi{t)ni{t,a — 0,x) = [2 — q{t,x)]v2{t)n2{t,a = l,x), 

^{t,a,x^x+{t))^0, ^{t,a,x^x+{t))=0, ^{t,a,x = x+{t)) = 0, 

ni{0,a,x) = ni{a,x), 71,2(0, a, x) = 712(0, x), n3(0,a,x)=0. 

(2) 
We have used the following notations 

• t the chronological time variable, 

• x: abscissa of transverse sections along the cephalic-caudal OT axis, < x < x+(t), 
X = cephalic tip, x = x+(t), caudal tip at time t, 

• < a < 1: degree of progression along the mitotic (M) phase or the intermitotic 
(Gi, S*, G2) phase, 

• Vi{t), (i = 1,2), homogeneous to an inverse of time, is the frequency of the corre- 
sponding phase duration, 

• D{x): "effective cellular diffusion" at section x, 

• < q{t,x) < 1: postmitotic neuronal population growth rate, or rate of differen- 
tiation (a more complete model should consider it as a nonlinear function of the cell 
densities) . 

The model should also be completed with a law for the OT elongation along the 
development, which means a description of the x+(i). It should naturally take the 
form 

-x+it)=F{niit),n2{t),n3{t)), (3) 



with a function F{-) that describes how the OT development is ruled by the different 
cell types. In this case the Neumann boundary conditions at x'^it) should be replaced 
by the no-flux conditions 

d 

x+(t)ni(t, a, a;_|_(t)) + D{x)-;r-ni(t, a,x — a;+(t)) — 0. 

ox 

To access such nonlinear coefficients as q{-) and F{-) is too demanding in a first 
stage, and thus simplifications are needed. Therefore we reduce the model complexity 
with two assumptions. Firstly, for days between E2 and E6 we can neglect asymmetric 
cell proliferation which leads to g = and 713 = 0. Secondly, we assume that the OT 
has constant length between days E2-E4 and EA-EG; this leads to fix x+ independent 
of time with a jump at the end of day EA. 

3 Reduced model for the symmetric proliferative 
phase without postmitotic migrating neurons {q = 0) 

We reduce the model complexity in the present work with two assumptions. Firstly, 
for days between E2 and E6 we can neglect asymmetric cell proliferation which leads 
to g = and n^ = 0. Secondly, we assume that the OT has constant length between 
days E2-EA and EA-EQ; this leads to fix x+ independent of time with a jump at the 
end of day EA. According to these assumptions, we consider only a simplified version 
of the model that we present now. We also present its discretization. 

3.1 The simplified model 

We define the region 

f]:={0<t<r, 0<a<l, 0<x<x+), 

with x+ constant and we also suppose that vi,V2, D constant. This leads to the model 
defined by the system of partial differential equations 



dni I dm _ n d rai 

rif 1" f^l r)n — J-^ 



dt ^ '''^ da ~ ^ dx^ 

8"2 I ,, 3"2 _ 7-| t?^Tt2 

dt ^ "^ da ~ ^~d^' 

completed with the boundary and initial conditions 
V2n2{t, a — 0,x) = vini{t, a = 1, a;), 
vini{t, a — 0,x) — 2v2n2{t, a = 1, x), 
^{t,a,x^ x+{t)) - 0, ^{t,a,x = x+{t)) = 0, 
ni(0, a, x) — n\{a, x), "-2(0, a, x) ~ ^2(0, x). 



(4) 



(5) 



For later purposes we impose the compatibility condition between the boundary 
and initial conditions: 

vin\{l,x) = V2nl{Q,x) (6) 

It implies that the solution is smooth enough for our analysis below, see [5]. 

3.2 Discretization by the method of finite diff'erences 

For every iiTGN, /gN, JgNwc define a mesh in Vl by the points 

{kAx,iAa,jAt), k = Q,l,...,K, i = 0,l, ...,/, j = 0,l,...,J 

with 

T . 1 . x+ 

At=—, Aa=-, Ax=^. 



For explicit schemes, these parameters have to be chosen with the stabihty (Courant- 
Friedrichs-Lewy) condition 

At At 

Aa (Ax)^ 

According to weh estabhshed methods [301 IS] , we approximate the derivatives in 
the differential equations by upwind (for first order derivatives) or centered (for the 
second order derivatives) differences on the mesh points. For a generic point of the 
mesh {kAx,iAa,jAt) we set as usual 

dni ^ ni{k + l,i,j) -ni{k,i,j) dni ^ ni{k,i,j) ~ni{k,i ~ l,j) 

dt ~ At ' da ~ Aa 

d'^ni ^ ni{k,i,j + 1) - 2ni{k,i,j) + ni{k,i,j - 1) / _ i r, 
dx^ " (Ax)2 ' " ' ■ 

In the following we will denote by: 

• n{k,i,j) value of the exact solution at the point (kAx,iAa,jAt) of the mesh, 

• u{k,i,j) value of the discretized solution at the point {k)Ax,iAa,jAt), 

and we consider the system of equations obtained upon application of these approxi- 
mations, with the corresponding boundary and initial conditions. 

The discretized model reads, for i = 1, ..., /, j = 1, ..., J — 1, k = 0, 1, ...,K, 

' u-i(k+l.i,j)—u-i(k,i,j) ^^ ui{k,i,j)—ui{k^i—lj) j-^ui{kA,j-\-l) — 2ui{k,i.j)-\-ui(k.i,j — l) 

~~Si l"^! AS ~ ^ ' (A^p ' ' 

U2{k+l,i,j)-U2{k,i,j) I U2{k,i,j)-U2(k,i-l,j} _ j~,U2(k,i,j+l)-2u2(k,i,j)+U2{k,i,j-l) 
At ^ ^2 A3^ - ^ ^A^p , 

V2U2ik,0,j) = ViUi{k,I,j), 

viui{k,0,j) = 2v2U2{k,I,j), 

. ui{0,i,j)^n^{i,j), U2{0,i,j) = n'^{i,j). 

(8) 
It is again convenient to also assume a compatibility condition between the bound- 
ary and initial conditions in the discretized model 

«in;(/,j)-«2n^(0,j) j = 0,l,...,J. (9) 

We use u to indicate both ui and U2, and v for vi and V2. Therefore, from equation 
([U, we deduce 

u{k + l,i,j)~u{k,i,j) u{k,i,j) -u{k,i~l,j) 

^ +V T = 

At Aa 

^ M(fc,i,j + 1) -2u{k,i,j) + u{k,i,j - 1) 

(Ax)2 
In other words, we obtain u{k + l,i,j) thanks to the iterations on the time label 

u{k + l,i,j) = u{k,i,j) ~v-—{u{k,i,j) ~u{k,i- 1J))+ 

At " (10) 

+-D7^-^("(fc,«,J + 1) ~2u{k,i,j) +u{k,i,j - 1)). 

We use the discrete Neumann boundary conditions 

u{k,i,0) — u{k,i,l), u{k,i,J+l)—u{k,i,J), (11) 

for fc = 0, l,...,i^, i = 0,1, ...,/. 

Applying the induction formula (|10p to calculate the values of ui and U2, and 
bearing (|lip and the initial conditions in mind, we obtain the values of the discretized 
solution 

ui{k,i,j), U2{k,i,j), fc = 0, l,...,ii:, 1 = 0,1,...,/, j = 0,l,...,J. (12) 

This describes the numerical method we use in the sequel together with experi- 
mental data. 



3.3 Convergence of the approximate solutions to the exact so- 
lution on the points of the mesh 

The convergence analysis of this numerical scheme is standard and can be found in 
[Sni IS]- We give here a fast account for the sake of completeness. 

We use the notation n{k,i,j) to indicate both ni{k,i,j) and n2{k,i,j), and also 
u{k,i,j) ior ui{k,ij) and U2ik,i,j). 

Applying Taylor's formula to the exact solution at the points of the mesh we obtain: 

n{k + l,i,j) -n{k,i,j) ri{k,i,j) -n{k,i- l,j) _ 

At Aa /, n\ 

n{k,t,j + l)-2nik,i,j) + nik,t,j-l) ^''^> 

{Axf + ' 

where for some constants A, B, C, we have 

\^<AAt + BAa + C{Axf. (14) 

From P^ we obtain 

n(fc + l,i,j) = n{k,i,j) -v^{n{k,i,j) -n{k,i~ l,j))+ 

+D^{n{k,i,j + 1) -2n{k,i,j) + n{k,i,j - 1)) + ^At. 

Next, we denote the error introduced by approximating n(k,i,j) by u{k,i,j) as 
w(k,i,j) = n{k,i,j) -u{k,ij) (16) 

From (fTU)) and ([T5)) . we obtain 

w(k + 1, i,j) = (l - t'fi - 2i?(^) Ak, i,j) + v^^w{K I - 1, i) 
+^(^("'(fc' «'i + 1) + <k, i,j - 1)) + *Ai. 
From the stability condition (jT]) and using (ITT)) , we obtain 

\w{k + l,t,j)\< (^l-v^^^2D^)\w{k,t,j)\+v^\w{k,i~l,j)\~^ 
+D^{\w{k,iJ + 1)1 + \w{k,t,j - 1)1) + |*|A< 

For every k we define 

Affc = max \w{i,j,k)\. 



(17) 



(18) 



«j 



From (IT51) . we obtain successively 



f At At \ At At , , , 

Af.,, < (l - ^^ - 2i? (^ j A/. + v-M, + 2i^^^A4 + 1*1 At, 

Mk+i < Affe + |*|At. 
Summing over fc = 0, 1, ..., m — 1, m < K 

771—1 m— 1 

X] ^^fc+1 < H A4+mAi|*|, 

fc=0 fc=0 

Afm < A^o +"iAt|*|. (19) 

Let us observe that u{0,iAa, jAx) = n{0,iAa,jAx), and therefore we have 

w(0, z, j) = w(0, ^Aa, jAa;) =0, (i = 0, 1, ..., / , j = 0, 1, ..., J). 



Bearing in mind that Mq — and At — ^, from (fTO)) we obtain 

T 
K 

Since ^ < 1, we also conclude that M,„ < |*|r. 

On the other hand, from (fT4)) . we finally obtain 

M„ < {AAt + BAa + C{Axf) T, m = 0, 1, ..., JiT. 
This proves that 

max\n{k,i,j) — u{k,i,j)\ ^ as Ai — >• 0, Aa — >■ 0, Ax ^ 0. 

The previous reasoning is valid for both components of the solution of the system ^, 
therefore 

max |ni(/i;, i, j) — ui(A:, i, j)| — > and max |ri2(fc, J, j) — U2(fc, i, j)| — > 0, 
as At ^0, Aa-^ 0, Ax -> 0. 

4 Dynamics of the simplified model from E2 to £"4 

We now study the dynamics of the model between E2 and EA. It is assumed that only 
symmetric mitosis takes place during that period. Then, we use the model (discretized 
by finite differences) given by equations ([5]). 

4.1 Experimental data recorded from the optic tectum at E2 

Records of number of cells counts performed on successive segments along the optic 
tectum cephalic-caudal axis were used as initial data for the model. Such informa- 
tion at day E2 corresponds to 16 segments, 25 nm. We have named these segments 
"sections". Every segment is identified by a value of the abscissa x. 

The initial model length is Le2 = 16 x 25/im = 400fj,m , and the total initial cells 
number is Ne2 = 800. 

The distribution of cells number per section is reported in TABLE 1 (Annexe) 

4.2 Determination of parameters 

Cell cycle duration between E2 and E4. According to experimental data, in the 
interval E2 — EA, the number of cells increases from approximately around 800 to 
around 9200 cells. Assuming that during this period the whole neuro-epithelial cells 
population exhibits a homogeneous proliferative behavior with a cell cycle duration 
around a mean value, then the averaged duration can be estimated as 13.6 hours. 

• Cell cycle duration: Gi -^ S + G2 + M — 13/i36mm = 13. 6/1. 

Intermitotic cycle duration Gi -t- S* -f G2 is 96 % the cell cycle duration; hence, 

• Duration of Gi -f S* -h G2 = 13.056ft. ^ 13ft03mm. 

The mitotic cycle duration M represents 4% the cell cycle duration : 

• Mitotic duration M = G.544/i ~ 33min. 

Number of cycles Gi + 5 + G2 + i\/ between E2 and E4. 

48feg ^0,0 

"^-^ ^ 13:6^1 = ^-^^ 

Total number of cells at E4. Using the formula 

Ne4 = Ne2 X 2"^2-4 



we obtain 

NEi = 800 X 2^-^^ = 9241 

Cellular density. It was observed experimentally that the average cellular density 
<^B2-4 between E2 and EA remained approximately constant, 3.85-^. Then, it will 
be assumed that between E2 and EA the model average cell density corresponds to 
this value. 

eel 



Model length at EA. 



<>E2-4 = 3.85- 



Ne4 9241ce/ „,„„ 

Le4 = 1—^ = r = 2400um 

Sp 3 85-££^ 

"£2-4 'J-OO^m 



Number of 25/im thickness sections comprising the model at EA. 

_ 2400^771 



400/xm 



==6 



Calculation of vi and V2 between E2 and EA. The following formulae define 
Inland V2, 



24hs 



• ^ — It Ui iilLUipilclOl.O t^i 1 <^ 1 1^2 P*^ 


^ interphase duration in /is 




24/is 


i;2 ■; J mitotic 


duration in hs 


In our case we have 




wi 13 


i"0.M4"^^-^2 


and we obtain 




vi ^0.54, 


W2 = 0.023. 



Time interval. We take t = at _E2 and t = T = A8hs at EA. We decompose the 
interval in sub- intervals of 10 minutes and we take At = lOrnm = ^hour. 
The index will be fc = 0, 1, ..., 48 x 6 = 288 , K = 288. 



Interval on the longitudinal axis. We decompose the interval [0,400] in 16 sub- 
intervals of length Aa; ~ 25fim. 
The index will be j = 0, 1, ..., 16 J = 16. 

Interval of variation of a. The interval of variation of a is [0, 1]. We will choose 
Aa to satisfy the stability condition 

At At At At 

^^A^ + ^^(A^^i' ^^a^ + 2^(a;^^^- 

Since 

1 1 (um)'^ 1 

vi = 0.54-, At = -h, L» = 6 ^^, , Ax = 25/im, V2 = 0.023- (20) 

h 6 h h 

the condition of stability will be the following one: 

0.54^ + 12l/^<l, 0.023^ + I2i/^<1. 

Aa 625 - ' Aa 625 " 

From the stability condition ([7]), it is sufficient to choose Aa such that 

1/6 1/6 

0.54-^ + 12-f- < 1, or Aa > 0.09. 
Aa 625 - 

We will adopt that index will bei = 0,l,...,10 that is 

7 = 10, Aa = 0.1. (21) 
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4.3 Induction formula for ui 

ui{k + l,i,j) = ui{k,i,j) - vi^{ui{k,i,j) - ui{k,i- 1, j)) + 

(22) 
+Dji^iuiik,i,j + 1) - 2ui(fc,i,j) + Ul(fc,^,J - 1)). 

Summing over i from, i = to i = / in (19), we obtain 

/ I (I I ^ 

X; ui{k + l,i,j) = Y. uiik,i,j) ~vi^< J2 uiik,i,j) - X) ■"i(^>«- I'j) ( + 

+^7^ Is wi(fc,»: J + 1) - 2 E ui(A:,z, j) + E ?/i(fc,i,j - 1)1 . 

U=0 i=0 i=0 J 

Observe that 

^ui(A:,i,j) -^ui(fc,i- l,j) ^ ui{k,I,j) -ui{k,-l,j), 
therefore, from (|23l) we obtain 



(23) 



X; ui(fc + l,i,j) == X) ui{k,i,j) -vi^{ui{k,I,j) - ui(fc, -1, j)} - 

+ ^7^ IE "l(fc:«>i + 1) - 2 E Ul{k,t,j) + Y: Ui{k,t,J - 1) 



(24) 



Li=0 i=0 i=0 

Let us also observe that the sum 

/ 

Ui{l,r) ^^ui{l,i,r) 

is the total number of interphasic cells (Gl + S* + G2) at instant / and at section r. 
Replacing in (|24| we obtain 

(k,j) -vi-^ |ui(fc,/,.7) - ui(fc, -1,.7)| + 

A* ..... .... . ... (25) 



f/i(fc + 1, j) = f/i(fc, j) - z;i|f {«i(fc,/, j) - «i(fc, -1, j)} + 

+DxBf {Uiik,j + 1) - 2 C/i(fc, j) + UiikJ ~ 1)} 



ui(fc, I,j) is the number of cells that complete the interphase Gl + S* + G2 at instant 
fc, section j . 

ui{k, 0, j) is the number of cells that begin the interphase Gi + S* + G2 at instant fc, 
section j. 

The number ui {k, 0, j) — ui{k, 0, j) depends on the total number of intephasic cells 
Gl + 5 + G2, at instant k and section j ; therefore, it is a function of Ui{k,j). 
It is difficult establishing this dependence law. It will be assumed that there exist a 
constant a(fc, j) , that depends on instant k and section j, such that 

ui{k,0,j)~uiik,0,j) = aik,j) Uiik,j). (26) 

Now we replace in (PSJ) , bearing in mind (PE)) and the values of f 1 , Ai , Aa; , Aa , D 
in (1201), to obtain 

C/i(fc + 1, j) = (1 + 0.9a(j))f/i(fc, j) + 0.0032 (cM^^iiil+CMM^ _ [/^(fc, j)) • (27) 

4.4 Approximation of the induction formulae 

Since, for each fc, the function Ui{k,j) is continuous in j, we conclude that 

0.0032 f^MM±l)±lMM^_^^(,^^.) 
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will be very small. Therefore, it will be ignored in a first approximation. Then, the 
formula (|27p takes the following approximate form 

C/i(fc + 1, j) = (1 + 0.9a{j)) C/i(fc, j). (28) 

Multiplying expressions I^E^ for fc = 0, 1, ..., 287 we obtain 

287 287 

nC/i(fc + l,j) = (l + 0.9a(j))288 nC^i(fc,j), or 

k=0 k=0 (29) 

Ui{A8hs,j) = (l + 0.9a(j))288c/i(0,j). 

Now, Ui{A8hs,j) represents the number of cells generated between E2 and E4: by 
the cells that at E2 were located at segment j. Given that during this period the 
model length increases 6 fold, it is assumed that, at Ei these newly generated cells 
will expand over six 25/im thickness segments. 

Thus, cells located at segment at E2 give origin to cells that at E4 occupy s - 
sections from s = 1 + 6{j — 1) to s = 6j . The values of j will be j = 1,2, ..., 16, 
corresponding to 16 sections of the model at E2. 

Ui{48hs, s) indicates the number of cells located at segment of abscissa s at EA. 
Ui{A8hs, sY^Zi.Qi -^i) will indicate the average number of cells in the segments of 
abscissa s at Ei, for 1 + 6{j — 1) < s < 6j. 
We have 

fT /AQh ^«=6J t/i(48fas,j) 
t^i(48/is,s).=i+60-i) = Q • 

For the sake of brevity, we will sometimes write Ui{4:8hs, s) instead of L^i(48/is, sY^~A_„, ._^^ 
for the corresponding value of j which appears in the above mentioned formula. We 
will have 



C7i(48/is,s) = i(l + 0.9a(i))288(7i(0,i), 1 + 6{j ~ 1) < s < 6j, or 
C7i (48/15, s) = mUi{0,j), J = 1, 2, ..., 16, 

where 

/3(j)^i(l + 0.9a(J))^«^ 
6 

We define the cell proliferation at segment j as the quotient 

Ui{A8hs,j) 
cel.prol.{j) = — -. 

Ui{Qhs,]) 

Since Ui{A8hs,j) = 6 C7i(48/is, s), we obtain 

cel.prol.{j) = 6/3(j). 

/3(j) is a coefficient proportional to the cell proliferation at segment j . 
The relationship between a{j) and I3{j) for j = 0, 1, ..., 16 is 



(30) 



<J) = ^ ([6/3(j-)]^ - l) m = ^ (1 + 0.9a(j))^ 



It will be assumed that proliferation is constant along the model, with a value equal to 
an average proliferation. This is equivalent to assume that a{j) — a and f3{j) — (3 
are independent of j. 

Finally, we will have 

C7i(48/is,s) = /3(j)t/i(0,j) J = 1,2, ...,16. (31) 
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We will calculate the constants a and /3 . We will consider a hypothetical segment j 
in which the number of cells will be obtained by calculating the average of the total 
number of cells of 16 sections in which the model was divided at E2. This number is 

^-- 

Likewise, we will consider that cells of section j originate a number of cells that at 
£'4 occupy 6 sections of the model. The number of cells at each one of the above 
mentioned sections is obtained by calculating the average of the totality of cells of the 
96 sections comprising the model. This average value is 

9241 „^ „^ 

= 96.26 

96 

and, therefore, 578 cells correspond to 6 sections. We will apply the formula (P^ 
taking 

t/i(0,j) = 50 f/^(48/is,i) = 578. 

We obtain 578 = (1 + 0.90)2^^50 

and therefore 

"578] 218 



. 50 

We conclude that 

a ^ 0.0095 and /3 ^ 1.93 (32) 

Replacing in (|5T]) we obtain 

t/i(48/is,s) = 1.93C/i(0,j) l + 6(j-l)<s<6j, i = l,2,...,16. (33) 

4.5 Number of cells given by formula (30) and comparison with 
experimental data 

The distribution of cells number per section is reported in TABLE 2 (Annexe) 

4.6 Approximation for linear interpolation of the number of 
cells in the intermediate sections and comparison with experi- 
mental data corresponding to £'4 

We will put 

C7i(s) = C7i(48/is, s), 1 + 6(j - 1) < s < 6j, j = 1, 2, ..., 16. 

We will use the formula of linear interpolation 

D 

forl + 6(j-l)<s<6j ,j = l,2,...,16. 

The value C/i(97) used in the calculation for the slope in interval 91 < s < 96 is 
obtained by means of the following trick. Let us observe that in the initial tip of the 
model at El the variation in the cell number between segments j — 1 and j = 2 is, in 
increasing sense 

C/i(0,j = 2)-t/i(0,j = l) = 7 

This induces us to extrapolate TABLE 2 by putting Ui{0,j = 17) = 42, after which 
the variation at the model final tip will be in diminishing sense, of the same magnitude, 

C/i(0,j = 17)-[/i(0,j-16) = -7 

Applying formula ((33| for j = 17 we obtain 

Ui{s) = Ui{A8hs,s) = 1.93 x 42 = 81 for 97 < s < 102. 

We take the value C7i(97) = 81. 

Total number of cells per section s at EA is reported in TABLE 3 (Annexe) and 

Fig.m 
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Figure 1: Total number of cells per section s in E4. 



5 Dynamics of the model between £"4 and E6 

5.1 Initial data 

The number of cells calculated by the model for an EA optic tectum was taken as 
initial data to calculate the evolution between EA and E6 and to estimate the number 
of cell each section will originate over that period. 

At EA the model is composed of 96 segments , 25/im thickness. 
The model initial length at E4 is Lea = 96 x 25/xm = 2400^to 
The total initial number of cells of the model at E4 is Ne4 = 9241. 

Table 3 illustrates spatially ordered distribution of number of cells/segment math- 
ematically (left) and empirically (right) estimated. 

5.2 Determination of parameters 

Cell cycle duration between EA and E6. Based on assumptions introduced in a 

previous section (4.2.) the average duration of the complete cell cycle G1 + S + G2 + M 

can be estimated as around as around 38 hs 15 min. 

Cell cycle duration d + S + G2 + M ^ 38.25ft. = 38ftl5mm. 

Assuming that the interphase Gi + S + G2 is 96 % the complete cycle, 

Interphase duration Gi + S + G2 = 36.72h = 36ft43mm. 

The mitotic length M is 4 % the complete cycle length, 
Mitotic duration M = 4% of 38.25/i = 1.53/i = lh32min . 

Number of cycles Gi + S + G2 + M between EA and E6. 

A8hs 



riEi 



= 1.255 



38.25/is 
Total number of cells of the model at E6. Using the formula 

Ne6 = Ne4 X 2"^4-6 

we obtain 

Nes = 9241 X 2^-255 ^ 22055 

Cellular density. It was experimentally observed that the average cell density SEi^e 
of the Tectum between EA and E6 remains approximately constant, with a value of 
3.85 eel/ iim. 
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We will suppose therefore that Se^^b — 3-85 eel/ ^m. 
Model length at ^6. We have Leg = #^^ = ■?lf'^?^ = 5729 rim 

Number of 25/im thickness segments comprising the model at E6. 

n"'' of segments =^t^ = 5729^™ ^ 230 

c> 2ofi7n 25/j.m 

Coefficient of cell diffusion between E4 and ^6. 

De.-. = ii^ = 2.387... - 2.Vm 

Calculation of ui and V2 between £'4 and £'6. The formulae that define vi and 
V2 are the following ones 

, J_ = „er of interphases d + ^ + G2 per day = .^t^rphase'duration ^n hs 

- ^ ~ n'^^ of mitosis M per day — ^^'"' 

We have — = ^^^ and — — -r$n therefore 

vi 36.72 V2 1.53 

„,.2|^.L53, .,.H5. 0.064. 

Time interval. We take t = at EA and t — T — AShs at £6. We decompose the 
interval in sub-intervals of 10 minutes each and we take 

At = lOmin = —hour. 
6 

The index wiU be fc = 0, 1, ..., 48 x 6 = 288, K = 288. 

Interval on the longitudinal Eixis. We decompose the interval [0, 2400] in 96 sub- 
intervals of length Ax = 25/irn. The index will be j — 0, 1, ..., 96 , J — 96. 

Interval of variation of a. The interval of variation for a is still [0, 1]. We choose 
Aa so that the stability condition ([7]) is satisfied 

At At At At 

iyi-— + 2D——-<l vi-— + 2D——-<l. 

Aa (Ax)2 - Aa (Aa;)^ " 

Since 

1 1 (um)'^ 1 

vi = 1.53-, At ^-h, D^ 2.4 ^^, ' , Ax = 25Aim, V2 = 0.064-. (35) 

h 6 h h 

the stability condition becomes 

:.53l^ + 4,8ig<l. 0.064^ +4.8lg<l, 
Aa 625 Aa 625 

and it is sufficient to choose Aa > 0.256. We adopt 

Aa = i (36) 

The index will be i = 0, 1, 2, 3 / = 3. 

5.3 Induction formula for Ui 

ui{k + l,i,j) = ui{k,i,j) -vi^{ui(k,i,j) -ui{k,i- l,j))+ 

(37) 
+DiEF)^{uiik,i,j + 1) -2ui{k,iJ) + ui{k,i,j - 1)). 

Proceeding along the same lines as in case of the model dynamics between E2 and 
£4, we obtain the formula 



Ui{k + l,j)= Uiik,j)-vi^{uiik,I,j)-u,ik,-lJ)} + 

+Dj^ {Ui{k,j + 1) - 2 Ui{k,j) + Uiik,j - 1)} , 
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where Ui{l,r) = X)i=o "i(^'*'^) ^^ ^^^ total number of interphasic cells Gi + S + G2 

at instant I and at section r, 

ui{k, I,j) is the number of cells at the end of the interphase Gi + S + G2 at instant 

k , and at section j, 

ui{k, 0, j) is the number of cells starting the interphase Gi + S + G2 at instant k , and 

at section j. 

As we did it for the model between E2 and EA, we suppose that there exists a 
constant a{k,j) such that 

ui(fc,0,j)-^i(fc,0,j)=a(fc,j)C/i(fc,j), (39) 

and that between EA and E6 a is also independent of fc, i.e. a{k,j) = <x{j). 

Replacing in ([55]) . and bearing in mind ([M)) . the values of vi, At, Ax, Aa, D in (|55|) 
and (|5S)) . we obtain 

C/i(fc + 1, j) = (1 + 0.765a{j)) Ui{k,j) + 0.00064 (^ c/i(fcj+i)+c/i(fc.j--i) _ Uiik,j)'^ . 

(40) 

5.4 Approximation of the induction formula 

The term 

fUi(k,j + l) + Ui(k,j -I) , / 

0.00064 f '^ '■' ' '^ '-^ '- - Ui{k,j) 

effectively acquires very small values which, in a first approximation, will be ignored. 

Formula (^0]) takes the following approximate form 

f/i(fc + 1, j) = (1 + 0.765a(i)) f/i(fc, j). (41) 

Multiplying expressions (|i(I)) for fc = 0, 1, ..., 287 we obtain 

[/i(48/is,j) = (1 + 0.765a(j))'''C/i(0,j). (42) 

Ui{4:8hs,j) represents the number of cells generated between E4 and E6 by the cells 
that at EA were located at section j . Given that during this period the model length 
undergoes a 2.4 fold increase, it is assumed that, at E6 these newly generated cells 
will expand over 2.4 segments of 25/iTO thickness. 

We will indicate with s the abscissa of transverse sections of the model at EQ. At 
each abscissa s = 1,2, ..., 230 will correspond a segment of 25/iTO thickness. 

It is assumed that the cells located at the j-th segment at EA give origin to cells 
that, at E6, occupy the part of the model corresponding to 2.4(j — 1) < s < 2.4 j, 
J = 1,2, ...,96. 

Ui {AShs, s) indicates the number of cells located at the segment of abscissa s at E6. 

f7i(48/is, s)^~2 4^ -i-) '^ill indicate the average number of cells at the segment of 
abscissa s at E6 , for 2.4(j — 1) < s < 2.4 j. We have 

Ui{A8hs,s)^^^/^^_^^ = — (43) 

The values of j will be 7 = 1, 2, ..., 96, which correspond to 96 segments of the model 
at EA. 
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In order to abbreviate notations, in the following paragraphs, the formula 
UiiA8hs,s) = UiiAShs, s)lZl%_r) 
will simply be written as [/i(48/is, s) for the corresponding value of j. 
We will have 

Ui{A8hs,s) = ^{l + Q.765a{j)f^^Ui{0,j) 
for 2.4(j - 1) < s < 2.4 j, j = 1, 2, ..., 96 or 

UiiA8hs,s)^mUiiO,j) 

where 

m^^a + 0.765a{j)f'' (44) 

We define the cell proliferation at segment j as the quotient 

Uii48hs,j) 
Ui[Qhs,j) 

Since 

Ui{A8hs,j) = 2.4 t7i(48/is,s) 

we obtain 

prol{j) = 2Am- 
The /3(j) coefficient is proportional to the cell proliferation at section j. 

Experimental data indicate that during the interval EA — EQ cell proliferation is 
not uniformly distributed along the cephalic-caudal axis. Thus, a coefficient j3{i) is 
introduced by means of the following simplified expression 

r 0,lj + 0.525 i/j = 1,2,3,4,5,6 

m={ 1-125 i/ J -6, 7, ...,48 (45) 

[ -0.01J + 1.605 i/ j ==48,...,95,96 

From (144p we obtain 

and replacing by (j45p we obtain the final expression for a{j). 
Finally, we will have 

[7i(48/is,s) = /3(j)C/i(0,j) for j = 1, 2, ..., 96, (46) 

where (i{j) is given by PSI) . 

5.5 Number of cells given by formula (43) and comparison with 
experimental data 

Number of cells per section s is reported in TABLE 4 (Annexe) 
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5.6 Approximation by linear interpolation of the cell number 
in intermediate segments of the model at E6 

We will denote by ?7i(48/is, s) = Ui{s) the number of cells at segment s of the model 
at E6 obtained by linear interpolation between values of TABLE 4. 

We will have 

Ui{s) = (l - '—^yF^) ^'^^-^^^ ^)) + '-^fi''h ,{2Aj) (47) 

for sections of the model in E6, with abscissa, 2.4(j — 1) < s < 2.4 j . 
We apply formula (|47| to calculate Ui{s) for integer values of s. 

Total number of cells per section s in EA is reported in TABLE 5 (Annexe) and 
Fig.m 
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Figure 2: Total number of cells per section s in E6. 



6 Discussion 

6.1 How to analyze the developmental role of cell proliferation 
in a developing system? 

The proliferation of NE cells is not only involved in generating the appropriate num- 
ber of cells in the CNS. A spatially organized cell proliferation activity also possesses 
morpho- and histogenetic effects [32j [23l |33j |34] . These developmental effects should 
depend on the temporal dynamics and the spatial organization of the proliferative 
activity. Taking this notion in mind, the present paper attempts at designing a math- 
ematical PDE-based model with the ability to describe the increase in cell number 
as a function of both time and space. The model emphasizes the occurrence of (a) 
changes in the cell cycle duration as function of time and (b) changes in the density 
of NE cells as a function of time and space. 

In our opinion this sort of PDE-based models can be useful tools to clearly un- 
derstand the cell proliferation spatial and temporal organization and to design formal 
models of cell proliferation that could account for its developmental effects. 

This paper is a first step towards designing a model that could help to globally de- 
scribe the temporal and spatial organization of the cell proliferation and to understand 
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its specific roles in the genesis of CNS patterning and its supracellular complexity. 

A current notion in theoretical biology is that the supracellular complexity is the 
net result of several DCBs operating simultaneously and interactively [2]. The cortico- 
genesis requires the integrated operation of several DCBs. Such organized operation 
requires spatially organized cell signaling processes which do not only control cell pro- 
liferation and migration but also establish cell-type differences and also determines 
tissue patterning i35j . It is currently accepted that these kinds of directive influences 
are mediated by gradients of specific signaling molecules, the so-called morphogens. 
A developmental gradient may be defined as an asymmetry in morphogen concentra- 
tion with a maximum at the site of the morphogen-secreting cells and decaying as a 
function of the distance from that site 3 . 

As an example, the non-uniform distribution of a morphogen in an embryonic field 
differentially determines the fate and phenotype of those cells [37] . In a same way, it 
can be proposed that the non-uniform spatial distribution of developmentally active 
signals infiuencing the proliferative activity should result in a non-uniform distribu- 
tion of the density of proliferating cells and also in the existence of space-dependent 
differences in cell density along a defined spatial axis. 



6.2 Stochastic approaches on the NE cell proliferation dynam- 
ics during the OT corticogenesis 

In the absence of directive influences a DCB should operate at random. In such condi- 
tion a biological time or space signal representing the DCB operation or representing 
a quantifiable direct developmental effect should display a random white noise-like 
dynamics. 

A previous work devoted to characterize the postmitotic neuronal migration dy- 
namics in the developing OT by means of standardized methods of frequency signal 
analyses [3S], proposes a mathematical expression that globally describes how cell 
migration operates in space. This expression includes deterministic and stochastic 
components. The deterministic component can be described as the sum of a linear 
global trend and a Fourier series while the stochastic component can be described as 
a non-correlated process (white noise). 

It is reasonable considering that, in the absence of directive infiuences the NE cells 
proliferation should operate at random and, consequently, the intercalation of newly 
generated cells into the whole system should be uniformly distributed. Under such 
condition the OT growth should proceed as a uniformly distributed tangential expan- 
sion; a characteristic that does not occur during the normal OT development. In a 
previous paper, aiming at determining whether the NE cell proliferation exhibits a 
spatial organization, spatial records of proliferating (mitotic) NE cells were analyzed 
within the framework provided by the stochastic point process model [H] . These anal- 
yses showed that the NE cell proliferation records can also be described as comprising 
two components: (a) a non-stationary one that can be deterministically described as a 
trend plus (b) a stationary non-correlated stochastic point process. The deterministic 
trend can be interpreted as manifestation of an asymmetrically distributed controlling 
infiuence that installs space-dependent differences in the NE cells proliferation rate 
along the cephalic-caudal axis. 

A more recent paper [39] indicates that the dynamics of a detrended NE cell pro- 
liferation record does not accurately coincide with a non-correlated stochastic point 
process with uniform probability distribution of intermitotic intervals. By contrast, 
stochastic point processes simulated with intermitotic intervals displaying lognormal, 
exponential and gamma distribution better approximate the NE cell proliferation 
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records than those simulated with a uniform probability distribution. 

6.3 The present conclusions on the temporal and spatial orga- 
nization of the NE cells proliferation during the OT corticoge- 
nesis. Morphogenetic correlations 

The PDE-based model presented in this paper was specifically design to obtain a bet- 
ter and more realistic characterization of the deterministic behavior exhibited by the 
proliferating NE cells. In our opinion, the present results appropriately describe the 
spatial organization of the proliferation activity and explain the differential growth the 
developing OT displays between E2E6. The model accurately estimates the increase 
in the total number of cells along the defined temporal windows {E2-EA and E4-E6). 
These increases closely coincide with values of cell density empirically estimated by 
means of cells counts performed on histological sections of developing OTs. 

With regards to the spatial organization of the NE cells proliferation, the model ap- 
propriately reproduces the empirically recorded position-dependent differences in cell 
density observed in successively ordered 25/i?7i length windows along the OT chepalic- 
caudal axis. In fact, figures 1 and 2 show that the cell density (number of cells / 
each 25fim windows) significantly changes as a function of their position along the 
cephalic-caudal axis. A comparison between these figures shows that, considering the 
OT as a whole, the NE cells proliferation significantly decreases as a function of the 
time (compare the profile of cell density values obtained at the end of the interval 
E2-EA with the profile of cell density at the end of the interval EA-EQ). 

It is interesting remarking that, apart from reproducing the increase in the to- 
tal number of cells, the profiles of cell density as a function of the space also allows 
accounting for some morphogenetic events. In fact, between E2E6 the OT under- 
goes a remarkable change in shape and size. The differences in cell proliferation as a 
function of their position along the cephalic-caudal axis should result in a differential 
intercalation of new cells within the OT wall along that spatial axis. The addition of 
fewer new cells at both extremes of the OT longitudinal axis should imply a lower OT 
growth rate at those sites. It must be noted that the low rate of NE cell proliferation 
at both ends of the cephalic-caudal axis reveals a low NE proliferation rate along the 
dorsal midline. This lower proliferation activity should also produce a lower growth 
rate along the entire dorsal midline. This assumption is corroborated by the fact that 
between E2E6 a deep medial groove appears along the dorsal midline. This change is 
a relevant morphogenetic effect since the appearance of a medial groove leads to the 
constitution of the right and the left OT hemispheres. 

All these data about (a) the changes in the number of cells as a function of the 
time, (b) the space-dependent differences in NE cell proliferation along the cephalic- 
caudal axis together with (c) the adequacy between the changes in size and shape 
that can be inferred from the model results, strongly suggest that the PDE-based 
model here presented can be appropriately and advantageously used to more precisely 
describe and to more deeply understand the developmental effects of a spatially and 
temporally organized NE cell proliferation during the OT corticogenesis. 
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